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' Abstract 

CO , Context. Methanol has a rich rotational spectrum providing a large number of transitions at sub-millimetre wavelengths from a 

range of energy levels in one single telescope setting, thus making it a good tracer of physical conditions in star-forming regions. 
Furthermore, it is formed exclusively on grain surfaces and is therefore a clean tracer of surface chemistry. 
J Aims. Determining the physical and chemical structure of low-mass, young stellar objects, in particular the abundance structure of 

Pm , CH3OH, to investigate where and how CH3OH forms and how it is eventually released back to the gas phase. 

^0 ' Mettiods. Observations of the Serpens Molecular Core have been performed at the lames Clerk Maxwell Telescope using the array 

receiver, Harp-B. Maps over a 4^5x5^4 region were made in a frequency window around 338 GHz, covering the Tk-^k transitions of 
methanol. Data are compared with physical models of each source based on existing sub-millimetre continuum data. 
Results. Methanol emission is extended over each source, following the column density of H2 but showing up also particularly 
^ . strongly around outflows. The rotational temperature is low, 15-20 K, and does not vary with position within each source. None of 

5— ( ' the Serpens Class sources show the high-A" lines seen in several other Class sources. The abundance is typically 10"' - 10"^ with 

^ ' respect to H2 in the outer envelope, whereas "jumps" by factors of up to IO--IO' inside the region where the dust temperature exceeds 

, 100 K are not excluded. A factor of up to ~ 10' enhancement is seen in outflow gas, consistent with previous studies. In one object, 

SMM4, the ice abundance has been measured to be ~ 3 x 10"'' with respect to H2 in the outer envelope, i.e., a factor of 10' larger 
than the gas-phase abundance. Comparison with C"*0 7=3-2 emission shows that strong CO depletion leads to a high gas-phase 
abundance of CH3OH not just for the Serpens sources, but also for a larger sample of deeply embedded protostars. 
^ Conclusions. The observations illustrate the large-scale, low-level desorption of CH3OH from dust grains, extending out to and 

beyond 7500 AU from each source, a scenario which is consistent with non-thermal (photo-)desorption from the ice. The observations 
also illustrate the usefulness of CH^OH as a tracer of energetic input in the form of outflows, where methanol is sputtered from the 
(N ■ grain surfaces. Finally, the observations provide further evidence of CHuOH formation through CO hydrogenation proceeding on 

grain surfaces in low-mass envelopes. 
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1. Introduction sitions observable in a single frequency window it is possible 

to obtain a coherent data set very efficiently, making methanol 

A long-standing goal of astrochemistry has been to determine a very suitable ti-acer of physical conditions, in particular- in 

5J] the physical and chemical conditions prevailing in star-forming low-mass star forming regions where the emission is optically 

Ci regions (e.g., van Dishoeck & Blake 1998). In this respect, dif- thin. Moreover, since the molecule is a slightly asymmetric top 

ferent molecules act as tracers of different physical components, mole cule, it traces very efficiently both density and temperature 

all depending on their formation history, their abundances, their (e.g., iMaret et aP l2005t lj0rgensen et al.l l2005bt iLeurini et aH 

chemical properties, etc. To effectively trace physical conditions 2007]). 

such as density and temperature over the large range of val- Methanol forms exclusively on ice-covered dust 

ues found m star-forming regions over the time-scale of star- ^^^^ ^^^^^^^ primarily through hydrogenation of CO 

formation, it is of great importance to have as many independent ( watanabe & Kouchi 2002; Fuchs et al. 2009). Observations of 

tracers as possible. Methanol (CH3OH), with its rich rotational interstellar ices show that methanol is indeed a prominent ice 

spectrum, is an excellent candidate tracing both temperature, component, with abundances of up to almost 30% with respect 

density, grain surface formation and energy injection simulta- sol id-state H,0 o r a fe w x 10"^ with respect to gas-phase H, 

neously during all phases of the early stages of stellar evolution. iPartois et all [19991: iGibb et al.1 iPontonmdan et all 

CH3OH is a slightly asymmetric top molecule with numer- 2004. In contrast, pure gas phase reac tions produce negl igible 

ous rotational transitions observable at millimetre- and sub- CH3OH abundances of less than 10"'" (iGarrod et al.ll2006l) . The 

millimetre wavelengths. Because of the large number of tran- question that naturally arises is how methanol desorbs from 
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the surface of a dust grain to be observed in the gas phase. Is 
it through thermal heating of the entire grain, or is it through 
non-thermal desorption, where cosmic rays, UV-photons or 
exothermic reactions provide local heating of the grain? The 
former mechanism is at play close to young stellar objects 
(YSOs), in the inner-most part of the molecular envelope where 
the gas temperature exceeds 100 K (van Dishoeck et al. 1995", 
Ceccarell i et aL 2000; van der Tak et al. 2000; Schoier et al. 
2002; M aret et al.l li2005; J0r gensen et al. 2005b) and in outflows 
where hot gas sputters the icy mantles (e.g. iBachiller et alJ 
119951; IBachiller & Perez Gutieireg Il997h . This mechanism 
allows for an efifective methanol enrichment of the environment 
and abundances are typically in the range of 10"^ to 10"^. The 
non-thermal mechanism dominates in co ld, dark clouds and 
the outer parts o f molecular envelopes ([H asegawa & Herbst 
1993HHerbst & C utJtjen 2006; Garrod etalJ lToO?; Oberg et al 



2009ah . Here, reported abundances typically have values of 
10"'° to 10"^. Hence, methanol also acts as a tracer of energetic 
processes in star-forming regions. 

So far most studies have concentrated on spectra at a sin- 
gle position or at most a few around them. Recently, large-scale 
mapping of weak molecular lines has become very efficient with 
the advent of array receivers such as the 16 pixel Harp-B re- 
ceiver on the James Clerk Maxwell Telescope (JCMT). This 
allows for a direct study of the entire protostellar system (pri- 
marily envelope and outflow) on scales of several arcminutes at 
15" resolution and for determination of density, temperature and 
energy input into the system. These observations will eventu- 
ally be compared directly to observations of another important 
grain surface product, H2O, to be done with the Herschel Space 
Observatory. By mapping the entire region, comparison with the 
different Herschel-beams (9"-40") will be straight-forward. 

The Serpens molecular core (also known as cluster A) 
is locat ed at a dist ance of 230±20 pc, following the discus- 
sion in lEiroa et aP feOOSi) . The Serpens molecular core con- 
sists of several deeply embedded sourc es, of which at least four 
are identified as contai ning protostars dWolf-Chase et al. I [19981 
iHogerheiide etaP ^^9). SMMl, SMM3, SMM4 and S68N. 
Large-scale continuum-emission studies have been performed 
to quantify the spectral energy distribution of all sources in 
order to classify their evolutionary stage as well as the dust 
prop erties of the molecular envelopes surrounding each source 
(e.g., 'Casali et al."1993'; Hurt & Barsonv 1996; 'Testi & Sargent^ 
1998; Davis et al. 1999; Larsson et al. 2000; Williams & Myers, 
20001) . most recently with the Spitzer Spac e Telescope as part 
of the Co res to Disks legacy program (c2d; iHarvev et al.ll200'7t 
lEvans et al. 2009). These studies show that three of the sources 
(SMM3, SMM4 and S68N) have relatively low bolometric lumi- 
nosities of ~ 5 Lq each, whereas SMMl has a higher luminosit y 
of ~ 30 Lq re.g.. iHogerheiide et al.l[T999t iLarsson et alJl2000h . 
The mass of each system (envelope and star) is in all cases less 
than 10 Mq. Recent interferometer observations by Choi (2009) 
show that SMMl is a binary system with a projected separa- 
tion of ~ 500 AU. The binary SMMlb appears less embedded 
than the primary (S MM la). This discovery has b een refuted by 
[Enoch et al.' ('2009') and Ivan Kempen et all (1200 9^ who both re- 
solve the disk around SMMl. Through detailed SED modelling, 
lEnoch et al.l (l2009t) finds a very high disk mass of ~ 1 Mq and 
that the inner parts of the envelope have been cleared out to dis- 
tances of500 AU. 

The region has been studied extensively at millime- 
tre (mm) and sub-millimetre (sub-mm) wavelengths in nu- 
merous molecular transitions (e.g ., iM cMullin et al. 1994; 

^ fT99i< 



ll999l:lMcMunin et alJl200(il) . however it was not included in the 
molecular surveys of isolated Class and I objects in Perseus 
and Ophiucus ( J0rgensen et al. 2004, 2005b; Maret et al. 2005*). 
The previous molecular studies conclude that SMMl, SMM3, 
SMM4 and S68N are all very similar to other young, low-mass 
stars with similar luminosities, such as IRAS 16293-2422 and 
NGC1333 IRAS4A and 4B in terms of abundances of sim- 
ple molecules that may be form ed directly in the gas phase, 
e.g., HCO+, CS, HCN CMcMulli n etan[T994l: iHogerheiide eTaH 
1 19991: iMcMulUn et all 12000 ). In Serpens, Httle has been done 
to quantify excitation and gas-phase abundances of molecules 
predominantly formed on dust grain surface s, like CH3OH, 
even though seve ral lines have been detected bv lMcMuUin et al] 
11991 l2000h and lHogerheiide et alJ 119991) . 

More direct observations of grain surface products have been 
made by Pontoppidan et al. (2004), who mapped infrared ab- 
sorption by molecules in the ice over a region extending 40" 
south of SMM4, but still located well within the molecular en- 
velope. The primary ice constituents were found to be H2O, CO 
(0.4-0.9 with respect to H2O), CH3OH (0.28 with respect to 
H2O) and CO2 (0.3-0.5 with respect to H2O; Pontoppidan et al] 
2008). The solid-state CH3OH abundance is one of the high- 
est reported to date, both when compared to H2O but also with 
respect to gas-phase H2 (3x10^^). At distances greater than 
12000 AU the CH30H-ice abundance drops beneath the detec- 
tion limit, corresponding to 3x10"^ with respect to H2. 

Besides the protostellar objects themselves, the region is 
permeated by large-scale outflows extending several arcmin- 
utes from the different sources with CO 7=2-1 velocities rang- 
ing from +10-15 km s~' with respect t o dt.sr - 8 - 8.5k m s~' 
dPavis etal.ll 19991 Graves et al. subm.). lGarav et alj (12002) ob- 
served the outflows from SMM4 and S68N in CH3OH 3k- 
2k emission and inferred CH3OH column densities of 1- 
2x10'^ cm"^, corresponding to molecular abundance enhance- 
ments of ~ 50-330, depending on outflo w position, consisten t 
with studies of other outflows (e.g.. .Bachiller et alJll995Lll998h . 

Here, the first map of rotationally excited methanol in the 
Serpens Molecular Core is presented of transitions which cover 
the energy range of E^p ~ 65 - 1 15 K. The paper is structured as 
follows. In Sect. |2]the observations are presented, and in Sect. |3] 
the observational results are provided along with radiative trans- 
fer modelling. Section|4]presents a discussion of the results, with 
a particular focus on the formation, desorption and excitation 
processes. Section|5]concludes the paper. 



2. Observations 

Observations of the Serpens molecular core were performed 
on June 20-22 2008 with the James Clerk Maxwell Telescop43 
(JCMT) on Mauna Kea, Hawaii. Observations were made of the 
1 k-6k rotational band of methanol {E^p ~ 65-115 K) at fre- 
quencies ranging from 338 to 339 GHz using the Harp-B ar- 
ray receiver consisting of 4x4 individual receivers dSmith et al.l 
,20031) . The telescope was pointed at two different locations near 
Serpens SMMl and two different locations near SMM4, see 
Fig. [Tj covering a total extent of 4!5x5!4. Observations were 
also made of C"^0, 7=3-2 at 329.330 GHz in a similar fash- 
ion, so as to be able to compare methanol and CO emission. 
Observations were made in jiggle-mode to achieve full spatial 



IWhite et al.1 119951: IWolf-Chase etaP fT998i: iHogerheiide et ah 



' The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research and the National Research Council of Canada. 
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Figure 1. Integrated emission, jTuBdv, of the strong 7o-6o A+ line at 338.41 GHz. Objects Serpens SMMl, 3, 4 and S68N are 
labelled. Offsets are in arcsec with respect to SMMl. Contours are at 3cr, 6cr, 9cr, etc. The units for colour are Kkms The 
pixel size is 7'.'5x7'.'5. The lines indicate outflow directions, with dashed representing blue-shifted outflow emission, and dotted 
lines representing red-shifted emission. The full white line extending south of SMM4 depicts the location of the four stars used by 
iPontoppidan et alj (12004 ) to measure the CH30H-ice abundance. Shown are also sample spectra at different positions, illustrated by 
triangles, from a single spatial pixel and provided in each inset in units of arcsec. The abscissa is for frequency in units of GHz and 
the ordinate is for main beam temperature in units of K. All spectra have the same scale to ease intercomparison. 



Nyquist-sampling using the Harp4 jiggle-pattern over a 2'x2' 
region. The beam-size of the JCMT at 338.5 GHz is ~ 15" and 
data were subsequently put on a map with a pixel-size of 7'.'5, 
i.e., one half beam-size to fulfill the Nyquist sampling criterion. 
The methanol observations were done using beam-switch with 
a throw of 180" in azimuth while the C"^0 observations were 
done using position switch to a clean position 1° off. The weather 
at the time of the observations was good with t(225 GHz) < 0. 1 . 
Pointing and calibration were checked at regular intervals, and 
the calibration error is estimated to be < 20% based on a compar- 
ison with standard calibration sources. Data were brought from 
the antenna temperature scale, T^, to the main beam tempera- 



ture scale, Tmb, by dividing with the main beam efficiency, ?;mb, 
which is 0.60 for the JCMT at these frequencies. The spectral 
resolution of the CH3OH spectra is 0.43 kms ', and the C'^O 
spectra were rebinned to the same value. 

The mean rms in the methanol spectra is 20 mK in 
0.43 kms"' velocity bins over the entire map. The Icr 
noise level on integrated emission has been determined as 
1 .2 (Trms \FWZIx6u wherc the factor 1 .2 accounts for the 20% 
telescope calibration uncertainty, cTi-ms is the rms noise, FWZI 
the estimated full width at zero intensity (taken to be 5 kms"') 
and 6v the velocity bin (0.43 km s"'). The 3cr uncertainty is typ- 
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ically ~ 0.1 Kkms ' per spatial pixel, but increases near the 
edges of the map by a factor of 2-3. 

Initial data reduction was done using the Starlink package. 
This consisted of putting the spectra on to a regular grid with 
pixel-size 7'.'5, subtracting linear baselines and co-adding all 
data-cubes. The Gildas package CLASS was used for subse- 
quent analysis. Standard data reduction revealed emission in the 
off-position for a subset of the methanol observations, but it 
was possible to remove this emission due to two facts: (1) The 
off-emission never coincided with on-emission, i.e., no double- 
peaked or "self-absorbed" line profiles were observed. In general 
the velocity offset between emission and absorption lines was of 
the order of ~ 5 kms"' compared to typical line-widths of 3- 
4 km s ' in the region where the off-emission was seen. (2) Data 
were recorded with a shared off-position, that is, even though 
the telescope jiggled on the on-position, no jiggling was done 
on the off-position. Therefore it was possible to recreate the off- 
emission in the strongest emission line at 338.409 GHz in the 
following way: first the strongest "absorption" hne was fitted 
with a Gaussian profile, then an artificial spectrum was created 
assuming a rotational temperature of 15 K, which was added to 
each spectrum. The spectra were then analysed and no absorp- 
tion features remained within the noise limit. 



3. Results 

In all spectra emission lines from the 1 k-^k rotational band 
of methanol are detected both on and off sources. In the cen- 
tral parts of each of the four YSOs up to eight emission lines 
are detected originating from ten transitions. The strongest line, 
the 7o-6o line at 338.409 GHz, has a peak brightness tem- 
perature of up to 0.5 K. All line profiles are Gaussian, no line 
asymmetries due to outflow activity or infall are observed. The 
FWHM of the emission lines are ~ 4-6.5 kms ' , consistent 
with other observations of me thanol line widths (e.g. lMaret et al.l 
120051; |j0rgensen et al.ll2005bh . Line widths change between ob- 
jects, but remain constant within each object, i.e., line widths 
do not change with position within an envelope. In Fig. [1] an 
overview is provided of the integrated emission in the strong 
7o-6() A^ line over the Serpens molecular core. The integration 
is numeric and has been done over the velocity interval -25 to 
-1-25 kms"' with respect to wisr = 8.0 kms"'. The four Class 
objects. Serpens SMMl, SMM3, SMM4 and S68N (labelled in 
Fig.lTJ are clearly seen. Representative spectra obtained at differ- 
ent positions sampling both YSOs, molecular envelopes and out- 
flows are shown. Besides the four Class objects, three distinct 
outflow knots are identified, SMMl-S, SMM4-W and SMM4- 
S. All of these bright knots coincid e with outflow po sitions as 
seen through CO observations (e.g., iDavis et af ] ll999h . Weaker 
structure is also seen in the map, for example an elongation of 
SMMl in the east-west directio n. This correspo nds to the direc- 
tion of a weaker CO-outflow (D avis et"anil999 ). The structure 
around S68N is more complex with the source itself being elon- 
gated in the NW-SE direction. Around this source a compact CS 
outflow has been disc overed at the same position angle as the 
elongation (Wolf-Chas e et al.lll998l) . Weak emission is also seen 
around S68N in the north-south and east-west directions. 

Representative spectra of the central parts of the four en- 
velopes are shown in Fig.|2]where emission is averaged over 3x3 
pixels (corresponding to 22'.'5x22'.'5 or 5600 AUx5600 AU). 
In the following, results will be presented first for the central 
parts of each envelope followed by results for the extended emis- 
sion in each envelope. Finally results will be presented for three 
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Figure 2. Spectra of Serpens SMMl, 3, 4 and S68N. 
The NGC13 33-IRAS2A spectrum is shown for comparison 
dMaret et"a D '2005). Spectra at the three outflow positions, 
SMMl-S, SMM4-W and SMM4-S, are also shown. Lines iden- 
tified in the CH3OH 1k-6k band are marked. Spectra were ob- 
tained by averaging emission from the peak pixel with its eight 
neighbouring pixels, i.e., in a box of 22'.'5x22'.'5. The SMMl, 3, 
4 and S68N spectra are shifted by +0.25, 0.5, 0.7 and 1.1 K, re- 
spectively, while SMMl-S, SMM4-W and SMM4-S are shifted 
by 2.4, 2.0 and 1.7 K. 



strong outflow knots in the map, named SMMl-S, SMM4-W 
and SMM4-S. 



3. 1 . Emission from ttie centrai part of tine enveiopes 

The emission lines arising from the central part of each enve- 
lope are fitted with Gaussian line profiles to obtain the integrated 
emission at higher accuracy. Results are tabulated in Table [1] 
along with the upper level energy, E^p, and line frequency. In 
the spectra presented in Fig. |2] the 3cr level is a factor of three 
lower, since they are each the average of nine spectra. Low-K 
lines up to = 3 are clearly detected, but high-ZT lines are not 
seen. This is in contrast to other low-mass YSOs with similar 
luminosities and/or distances and evolutionary stages, such as 
IRAS 162 93-2422 (e.g.,van Dishoeck et al...l995,) and NGC1333 
IRAS2A dMaret et al.li2005l) ! 

The integrated intensities are used to make rotational dia- 
grams for each of the four objects. These are shown in Fig. 
|3] Within the uncertainty, the logarithm of the derived upper 
level column densities, A^, divided by the statistical weights, g, 
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Table 1. Integrated line intensities, J Tub dy (Kkms"') over the central 22'.' 5 x22'.'5 of each source" and emission-weighted, 
average line-widths. 



Transition 


V (GHz)" 


Eup (K)" 


SMMl 


SMM3 


SMM4 


S68N 


SMMl-S 


SMM4-W 


SMM4-S 


+0E 


338.1245 


78.1 


0.25 


0.22 


0.30 


0.85 




0.09 


041 


-IE 


338.3446 


70.5 


0.58 


0.45 


1.09 


2.90 


0.52 


0.90 


1.38 


-hOA 


338.4087 


65.0 


0.71 


0.62 


1.35 


3.66 


1.06 


1.26 


2.21 


-2A 


338.5129 


102.7 


0.06 


0.01 




0.24 


0.08 


0.20 




±3A'- 


338.5420 


114.8 


0.11 


0.16 


0.15 


0.40 








-HlE 


338.6150 


86.1 


0.16 


0.08 


0.18 


0.45 




0.34 




+2A 


338.6399 


102.7 


0.09 


0.07 




0.26 








±2E' 


338.7223 


89.1 


0.21 


0.16 


0.28 


1.05 




0.13 


0.11 


Av (kms"') 






3.4±0.9 


4.4±2.5 


4.0+0.7 


7.2±1.3 


4.3+0.9 


4.0+0.3 


4.1+0.3 



" Non-detections are marked by ". . . ", for which an upper limit of 0.03 Kkm s"' applies. 

* Frequencies and upper level energies are from the Cologne Database for Molecular Spectroscopy ("CDMS; lMiiller et alj|2001h . 
Lines are blended. The frequencies and upper level energies are average values. 
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Figure 3. Rotational diagrams with data obtained from the four 
spectra shown in Fig. |2] The best-fit straight line is shown in 
each case as a dotted line. Error bars are for Icr uncertainties. 



fall on a straight line when plotted versus the upper level en- 
ergy divided by the Boltzmann constant, ks- No difference be- 
tween A- and E-type methanol is found, thus the abundances are 



equal, ;t;(E-CH30H) = .x:(A-CH30H). In the cold-temperature 
limit (T ~ 10 K) the abundance ratio is expected to he E/A = 
0.69 whereas it is 1 at higher temperatures, and thus littl e or n o 
difference in abundance is expected (e.g., Friberg et al. 19881) . 
The derived rotational temperature for each of the four envelopes 
is low, 15-20 K, significantly lower than the rotational temper- 
ature of ~ 80 K inferred for IRAS 16293-2422 and NGC 1333 
IRAS2A (Ivan Dishoeck et aI.lll995HMaret et al.ll2005 l). The esti- 
mated column densities are all of the order of ~ lO'^-lO'^^ cm"^. 
Results are tabulated in Table |2] Rotational tempe ratures and 
colum n de nsities are simi l ar to t hose found by iMaret et al.l 
(I2005h and lJ0rgensen et al] ( l2005bl) for a larger sample of low- 
mass. Class objects, with the exception of the two sources men- 
tioned above. 

3.2. Extended emission from the envelopes 

To quantify the spatial distribution of methanol in each molec- 
ular envelope, the radial distribution of emission from the five 
strongest lines is plotted in Fig. |6] Methanol emission is ex- 
tended in all envelopes, and the FWHM of emission is around 
2-3 times the 15" beam. Extended methanol has previously 
been reported in the en velope of the isolated Class object 
L483 dTafalla et al.ll2000 ) in the 2o-l() A+ and 2_i-l_i E lines 
(iiup - 12-20 K), where it was found that the methanol emission 
traces total gas column density. To verify whether this is true for 
the Serpens envelopes, it is necessary to quantify the physical 
structure of the envelopes. This is done through a combination 
of dust continuum emission at 850 fim as observed with SCUBA 
on the JCMT and modelling of the physical structure. 

To obtain the absolute abundances, three different methods 
are used. First, the CH3OH column density obtained from ro- 
tational diagrams is compared to the column density obtained 
from dust emission at 850 jim in the same beam. Second, the 
CH3OH column density from rotational diagrams is compared 
to the column density predicted from a physical model of each 
source. Third, a radiative transfer model is compared directly to 
the observed spectrum to constrain the abundance. Results are 
Usted in Tableg] 

3.2.1. Dust continuum emission 

It is possible to directly estimate the average column density over 
the central area of each envelope by assuming that the dust emis- 
sion is optically thin, thermal emission at a single dust tempera- 
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Table 2. Rotational diagram results for the observed methanol 
emission. 



Source 


r.-o, (K) 


N (10"* cm-^)" 


SMMl 


18.9+0.9 


4.1±0.7 


SMM3 


19.5±1.1 


3.4+0.7 


SMM4 


15.1±0.4 


12+1.4 


S68N 


14.4±0.2 


34+1.7 


SMMl-S 

SMM4-W 

SMM4-S 


12.6±2.4 
16.5±1.3 
9.3±1.0 


'-5.7 

4.4+1.4 
34+27 



" Column density is the sum of A- and E-type methanol column 
densities. 



ture, T^: 



N = 



Syd^ 



= 4.90 X 10"" cm-2 Sy{Jy)d(pcf 



/ 16.9K , 



(1) 



Here Sy is the continuum flux integrated over 22'.'5x22'.'5 at 
850 ;um, d the distance to the object, Ky the dust opacity at 
850 fim taken from Ossenkopf & Henning ( 1994), By the Planck 
function at a dust temperature and A the total surface area, 
in this case 7.1x10-'"' cm^. Finally a gas/dust ratio of 100 is as- 
sumed. 

Reduced data from the SCUBA Legacy archive have been 
used (Di Francesco et al. 2008) to estimate the integrated con- 
tinuum flux at 850 yum. An average dust temperature of 20 K was 
used following global estimates of the dust temperature over th e 
entire Serpens molecular core of 20 K from lSchnee etaP (l2005h . 
If the temperature changes to 10 K, the column density estimate 
is increased by a factor of 3.3, while if it is 30 K, the column 
density is decreased by a factor of 1.8. 

By using a dust temperature of 20 K, the average column 
density over a 22'.' 5 x 22'.' 5 region is in the range of 1^ x 
10^^ cm"^. This leads to CH3OH fractional abundances in the 
range of 1-20 x 10"^ with respect to H2, see below. 

3.2.2. Dusty modelling of envelope properties 

Since the dust (and gas) temperatures change near the YSOs, the 
CH3OH emission can be further quantified through modelling of 
the physical parame ters of the envelopes. This i s done following 



tne pnysicai parame ters 01 tne envelopes, inisjs aone roiiowmg 
ISchoier et alj ('2002') and J0rgen sen et al] (2002), where the dust 
continuum emission is modelled with a spherically symmetric 
envelope with a power-law density structure he ated from the in- 
side b y the protostar, using the ID code Dusty (Ivezic & Elitzuil 
Il997h . Dust opacities tabulated by Ossenkopf & Henning ( 1994) 
are used for densities of 10^ cm"^ and thin ice mantles (corre- 
sponding to the values in their Table 1, column 5). The central 
heating source is taken to be a black-body radiating at a temper- 
ature of 5000 K, but model results are not very sensitive to this 
parameter The output from Dusty is dimensionless, so to make 
an absolute calibration it is necessary to do so against the abso- 
lute luminosity and a given distance of the source. The output 
can be in the form of a spectral energy distribution (SED) and 
a radial profile of continuum emission at a user-specified wave- 
length. 

Only dust heating at temperatures below 250 K is considered 
(corresponding to a peak wavelength of ~ 12 //m). Thus there is 



a "hole" in the inner part of the envelope extending out to a ra- 
dius, Tin, where no attempt is made to model the dust emission. 
The physical extent of the envelope is defined by the parameter 
Y = nn/'"oiit- The density profile of the envelope is described by a 
power-law such tha t w(r) oc r^P. Three obs ervational constraints 
are used following |j0rgensen et al.l (|2002|) to identify the best- 
fitting dust model: the SED and the spatial distribution of con- 
tinuum emission as recorded by SCUBA on the JCMT at 450 i-im 
and 850 ^m, both availab le through the SCUBA legacy archive 
(iDi Francesco et al.ll2008h . 

SEDs were assembled from the literature, using data 
points from MIPS on the Spitzer Space Telescope, ISO- 
LWS and IRAS along with ground-based single-dish con- 
tinuum measurements at sub-mm, mm and cm wave- 
lengths. For this work, the SED database available at 
[http : //astro . kent . ac . uk/protostars/ proved invalu- 
able (IFroebrichI 120051) ! A full table of SED points is included 
in Appendix A. 

The reduced SCUBA continuum emission maps from the 
SCUBA legacy archive have been smoothed to equivalent 
beam sizes of 11" and 19'.'5 at 450 and 850 fim, respectively 
(iDi Francesco et al.l 120081) . Along with emission maps, maps 
containing the error on each pixel are also provided, facilitat- 
ing error analysis. The radial profile of each object was com- 
piled by considering the emission in annuli extending from the 
pixel containing the maximum value of emission. The maps 
were carefully checked, and directions containing contamination 
from other objects were avoided. 

Three Dusty input values were varied: the physical size of 
the envelope, Y, the opacity at 100 fim, r(100 fim) and the power- 
law slope of the radial density profile, p. A small grid of simu- 
lations was calculated and the best-fit model for each envelop e 
was found using a^^ method, following |j0rgensen et al.l (l2002h . 
The best-fit model parameters are listed in Table [3] along with 
physical parameters such as the extent of the envelope, the local 
density at the inner radius, the column density and the total en- 
velope mass. The best-fit model predictions are overplotted on 
actual data in Fig.|4] 

In general agreement between the modelled and observed 
properties are good, except in the case of the modelled 850 jum 
radial profile of S68N. The model prediction is consistently a 
factor of 2-3 lower than observed values. Part of the reason is 
that since several sub-mm sources are present close to S68N, 
only emission from the NW direction is modelled here, imply- 
ing that the uncertainty on individual data points is higher than 
in the case of SMMl, SMM3 and SMM4. In this particular di- 
rection there may still be contamination, and thus emphasis is 
placed primarily on modelling of the SED. Moreover, the ra- 
dial emission profile at 450 yum is much steeper than the 850 //m 
profile, indicating that there may be more cold, ambient cloud 
material toward S68N. The radial profile obtained at 450 yum is 
well reproduced by the best-fit Dusty model. 

From the Dusty model it is possible to estimate the average 
column density over a region of 22'.'5 x 22'.'5 and compare this 
to the CH3OH column density. The average column density is 
found to be ~ (1.0 ± 0.2) x lO^"* cm"^. This leads to an abun- 
dance of ~ 3 X 10"' with respect to H2 for SMMl and SMM3 to 
~ 30 X 10-'' for S68N (see Table|4]for details). 

3.2.3. Ratran simulation 

In the following the abundance is constrained through direct ra- 
diative t ransfer modelUng. This has been done using the code, 
Ratran dHogerheiide & van der Ta]jl2000l) in conjunction with 
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Table 3. Best fit Dusty model input parameters and results for the physical structure of the four molecular envelopes in Serpens. 



Property 


SMMl 


SMM3 


SMM4 


S68N 


IRAS2A 




37 


2.0 


6.3 


3.6 


21 


rf(pc) 




230 






250 


P 


2.1 


1.6 


1.8 


1.8 


1.9 


t(100 yum) 


3.8 


1.1 


3.0 


3.4 


1.8 


Y (=n„/ro„,) 


400 


1000 


1000 


800 


500 


'•in (AU) 


49 


14 


17 


18 


31.3 


ro„, (10* AU) 


2.0 


1.4 


1.7 


1.4 


1.6 


«,„ (IC cm-3) 


1.7 


0.94 


2.8 


3.0 


1.0 


"1000 AU (10* cm-3) 


3.0 


1.0 


1.9 


2.2 


1.4 


MHz) (lO^"* cm-2) 


1.1 


0.33 


0.89 


1.0 


0.53 


Envelope mass, T > 10 K (Mq) 


4.1 


2.5 


3.9 


3.6 


2.2 



-14 










- 16 










-18 




Y = 


400 




-20 




P = 


2.1 




-22 


ySMMl \ 




= 3.8 




-14 










^-16 








\* 


B 




Y = 


1000 




" -20 




P = 


1.6 




-22 


/ SMM3 \x 




= 1.1 




^.-14 










^ -16 










^-18 




Y = 


1000 




t!lO 






1.8 




B "20 




- P = 




-22 


. /sMM4 \ 


"^100 


= 3.0 




-14 










-16 










-18 




Y = 


800 




-20 


/ ^\ 


- P = 


1.8 




-22 


. /S68N \ 




= 3.4 
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0.100 
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Figure 4. Dusty modelling results of the molecular envelopes of 
SMMl, 3, 4 and S68N. The left panel shows the SEDs of all four 
objects, whereas the right panel shows the no rmalized radial dust 
emission profile as obtained with SCUBA (iDi Francesco et all 
120081) . Best-fit results are over-plotted as solid lines, and the val- 
ues for the best-fit models are given in the right-hand panel. 



molecular data from LAMDA (iSchoier et alj|2005h . This op- 
portunity was used to update coUisional data in the LAMDA 
database for methanol using newly released rate coefficients for 
collisions between para-H2 and A- and E-type methanol up to 
and includin g J — 9, corresp onding to an upper level energy of 
~ 500 K (Pottage et al.ll2004 ). The rate coefficients are given for 
temperatures from 5 to 200 K. The coefficients are not complete 
up to upper level energies of 500 K since they are limited to 
J < 10. In the case of A-type methanol with K - Q, this cor- 
responds to an upper level energy of ~ 100 K. Data have been 
extrapolated so that they are complete up to an upper level en- 



ergy of 385 K corresponding to the first torsional state. For the 
extrapolation only transitions between /^-ladders and A/ = 1-5 
are considered. The coUisional rate coefficients are found to 
be proportional to V?^ /A/ which is used for the extrapolation 
(Leurini et al. 2007). Level energies and Einstein A-values are 
taken from t he Cologne Da tabase for Molecular Spectroscopy 
(CDMS; MulleretalJl2"00lh . 

The density and temperature profile of the Dusty modelling 
is used for the physical structure of the envelope. No effort is 
made to model the line widths, but instead a constant, turbulent 
line width of 3.0-5.0 kms ' is assumed, depending on source. 
This is not a measure of the actual turbulence in the region, but 
only a means to ensure that the modelled linewidth corresponds 
to the observed. For each source a small grid of Ratran models 
has been run with a so-called "jump" abundance structure. The 
jump was located at rdu.st =80 K (7? = 50-125 AU) with an in- 
crease in the abundance of a factor of 1-10^ corresponding to 
most CH3OH evaporating off of the grains. When running the 
Ratran models, great care has been taken to ensure that the cen- 
tral region containing the jump is properly sampled. In order to 
do this the pixel-size is set to 0'.'2 (45 AU) and the central pixels 
were oversampled by a factor of 50. 

A method is then used to determine the best-fit model 
integrated emission compared to observed integrated emission. 
This is done separately for A- and E-type CH3OH to examine 
relative abundances. No significant difference is found between 
A- and E-type CH3OH abundances, thus confirming the results 
obtained from the rotational diagrams. All abundances are listed 
in Table |4] Diagrams illustrating the variation of with inner 
and outer abundance are shown in Fig. |5] Here it may be seen 
that the inner abundance is not well constrained, and it is only 
possible to provide upper limits ranging from 2x10"'' (SMMl) 
to 5x10"^ (SMM3 and S68N). On the other hand, the outer enve- 
lope abundance is very well constrained for all sources and lies 
in the range of 10"'' (SMMl) to greater than 10"** (S68N). The 
methanol enhancements, Xin/xout, are from < 2-3 (SMMl) to < 
200 (SMM3). Due to the lack of h\gh-K lines it is not possible 
to further constrain the inner abundance. 



3.2.4. Comparison of abundance measurements 

The three methods for determining the abundance may be di- 
vided into the following categories; 



1. Direct observational abundance assuming LTE in terms of 
CH3OH excitation and optically thin dust emission. 



SMMl 
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SMM3 SMM4 S68N 




-9 -8 -7 -9 



-7 -9 -8 -7 -9 
log[x(CH OH), J 



Figures. Reduced distribution as a func- 
tion of inner and outer abundance for A- and 
E-type CH3OH as well as the combined value. 
Contours are for 1, 2, 3, 4 and 5cr. The dot- 
ted line corresponds to the abundance derived 
using Niy^T) from Dusty modelling, and the 
dashed line is the abundance derived using 
A'(H2) directly from SCUBA dust continuum 
emission. The minimum x''^^^ is given in each 
panel. 



Table 4. Average total gas column densities and fractional abundances of CH3OH. 



Source 


"dusty 
(10" cm-2) 


"scuba 
(10" cm-2) 


X(CH3 0H)pyj.j,y 

xlO-^ 


xlO-^ 


xlO-' 


.^(CHsOH);;^' 

xlO"** 


SMMl 


1.3 


3.5±0.3 


3.6±0.9 


1.3±0.4 


1 


<0.3 


SMM3 


0.83 


2.0+0.13 


2.8±0.7 


1.2±0.3 


4 


<50 


SMM4 


1.1 


1.2±0.07 


13±3 


12±3 


5 


<5 


S68N 


1.2 


1.7±0.10 


28±2 


21±2 


10 


<40 



" Gas column density, ^(Ht), obtained from Dusty modelling over a region of 22'.'5x22'.'5 assuming a gas:dust ratio of 100. 
' Gas column density , A'(H2), obtained from SCUBA emission at 850 yum assuming a constant dust temperature of 20 K. 



^ Fractional CH3OH abundance summed over A- and E-type CH3OH. 
Inner and outer abundance in the Ratran jump-model. 



2. Combined direct and model abundance, still assuming LTE 
in terms of CH3OH excitation but modelling the dust emis- 
sion with a range of temperatures. 

3. Modelling of the CH3OH emission within a physical, spher- 
ically symmetric model with varying dust and gas tempera- 
tures and densities. 

Despite the obvious differences and levels of complexity be- 
tween the three methods, the results in terms of average 
methanol abundances are remarkably similar. Typically, varia- 
tions are within a factor of 2-3. 

In the following only results from the Ratran modelling 
will be discussed. Since this method takes the full density and 
temperature structure into account it will be the more accu- 
rate of the three as demonstrated for other molecules (e.g . , 
iHogerheiide & van der Tal3l2OO0l: |j0rgensen et ani2002l 12004^ . 

3.3. Outer envelope abundance structure 

In Fig. |6]the column densities as predicted by the Dusty mod- 
elling and the 850 //m dust continuum emission radial profiles 
are overlaid. The methanol emission profiles follow both the dust 
emission and column density profiles very closely over scales 
from ~ 5500 AU to greater than 12 000 AU. Because the distribu- 
tion is independent of rotational line, the rotational temperature 
(corresponding to line ratio) is constant throughout the envelope 
and the methanol column density scales directly with the enve- 



lope column density, similar to the case of L483 dTafalla et all 
2000!). Thus, over this range of radial distances, the average 
CH3OH gas abundance is constant. In the case of S68N the 
agreement between the observed and simulated radial profiles 
is not good since the simulation under-estimates the emission at 
larger radii. Furthermore, the outer envelope abundance as de- 
rived from rotational diagrams or from comparison with dust 
emission agrees to within a factor of 2 with the outer abundance 
as derived from the Ratran modelling. 

3.4. Emission from outflow positions 

Methanol is detected along the outflows throughout the emis- 
sion map (Fig. [TJ. In the following, focus will be placed on 
emission from three distinct knots labelled SMMl-S, SMM4- 
W and SMM4-S on Fig. [T] The elongated emission associated 
with S68N indicates that a compact shock is responsible for 
sputtering CH3OH from the grain mantles into the gas phase. 
This outflow has already been observed in CS 7=3-2 emission 
(Wolf-Chase et al. 1998). The line profiles associated with this 
outflow are symmetric and Gaussian in nature making an accu- 
rate disentanglement of the contributions from envelope and out- 
flow difficult. Unlike the above mentioned outflow knots, emis- 
sion is spatially coincident with the source position. For this rea- 
son the S68N outflow is not included in the following discussion. 

To investigate the properties of each outflow knot, spectra are 
extracted from the map and rotational diagrams are made in the 
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SMMl SMM3 SMM4 S68N 
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Figure 6. Spatial distribution of methanol emission in the four objects, Serpens SMMl, 3, 4 and S68N. The integrated emission 
of the five lines 7-6 +0E, -IE, +0A, +1E and ±2E is shown (black). The error-bars are for Icr uncertainties. Overlaid is the dust 
continuum emission as obtained from SCUBA (red) and the Dusty dust continuum modeling column density (blue). All profiles 
have been normalized to 1 . The normalization factor for the peak emission is shown in each diagram. The 3cr noise level is shown 
as a dashed line, and the grey, dashed line shows the beam profile. 
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Figure?. Outflow emission in Serpens as traced by the CH3OH 
7o-6() line. Contours are for red- and blue-shifted outflow 
emission and are at 3cr, 6cr, 9cr, . . . Crosses mark the positions 
of SMMl, SMM3, SMM4 and S68N. 



same manner as above (see Figures|2]and|3]l. The intensity of the 
strong 7o-6o line is comparable to the emission from SMM4 
in all three knots. The detected transitions and associated intensi- 
ties are listed in Table[T] The lines are shifted by ~ 1 km s ' with 
respect to each other and with respect to the source velocities. 
Line profiles are symmetric and Gaussian in nature, with no ob- 
vious line-wings. Nevertheless, methanol clearly traces the out- 
flow activity associated with SMMl, SMM4 and S68N as seen 
in Fig. |7] whereas little emission is associated with SMM3. In 
Fig.|7]the contributions from red- and blue-shifted emission are 
estimated by integrating over the velocity intervals from -5 to 
7 kms"' and 10 to 20 kms"', respectively. This is to be com- 
pared with the average source velocity, wisr = 8.5 kms"'. 

Results obtained from the rotational diagrams are tabulated 
in Table 12] The rotational temperature is ~ 10-15 K, indicating 
sub-thermal excitation at all three outflow positions, as found 
typically also in other outflows because of the relatively low 
dens ity compared with the central envelope (e.g.. lBachiller et aU 
The column densities range from a few X 10^'* cm ^ 
(SMMl-S and SMM4-W) to ~ 3x10'^ cm^^ (SMM4-S). To ob- 
tain accurate abundance estimates of methanol in shocked re- 
gions, it is again imperative to estimate the gas column density 
accurately. In this case '^CO 7=3-2 data from the JCMT science 
archive, obtained in the context of the JCMT Gould Belt legacy 
survey, are used (Graves et al. subm.). The principle for obtain- 
ing the column density from '^CO in the shocked gas is that 
emission from the higher velocity line-wings is optically thin. 
This assumption has been s hown to be valid for a number of 
molecular outflows (see e.g.. lBachiller & Tafallal[T999L for a re- 




Velocity [km s '] 



Figures. Comparison of outflow line profiles. For each out- 
flow position, SMMl-S, SMM4-W and SMM4-S, line profiles 
of i^CO 7=3-2 (red), C'^O 7=3-2 (blue) and CH3OH 7o-6o A+ 
(black) are shown. C'^O and CH3OH spectra have been multi- 
plied by 3 and 25, respectively, for clarity. SMM4-S and SMM4- 
W spectra have been shifted by 8 and 20 K. The masks at the 
bottom indicate the velocity range over which '-CO emission is 
integrated to obtain the outflow column density. 



view). In Fig. [8] the line profiles are shown for the three outflow 
positions analysed here. Emission from SMMl-S and SMM4- 
W is blue-shifted, while emission from SMM4-S is red-shifted. 
The CO line-wings are not as pronounced as in other well-known 
outflow sources (e.g., Blake et al. 1995), but are slightly shifted 
(up to +5 kms"') with respect to a yisr of 8.5 kms"'. This in- 
dicates that the flow s may be caused by J-type shocks rather 
than C-type shocks (iHollenbachI 1 1 997b . Molecular emission is 
expected to have a Gaussian velocity profile around the shock 
velocity in a J-type shock, whereas emission caused by C-type 
shocks is expected to show a "classic", triangular line profile. 
However, since the observed lines are only shifted by +5 km s ' 
at most, so the J-type shock speed would be lower than this. 
CH3OH is also efficiently destroyed in the sputtering process in 
J-type shocks at velocities greater than 10 kms '. The actual 
CH3OH enhancement may be lower in the outflows in Serpens 
compared to e.g., LI 1 57, somethirig whi ch is also discussed for 
the case of NGC2071 dGarav et al.l(2000h . Another possibiHty is 
that the flows are C-type shocks moving very close to the plane 
of the sky, in which case the line-wings will not be prominent. 
This is suported by SiO 7 = 2-1 observations of SMM4, where 
the lines were also found to be narrow (FWHM ~5 kms"'; 
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iGarav et"ani2002h . With the present data-set it is not possible 
to differentiate the two scenarios. 

Emissio n was integrated over the velocity intervals given in 
iDavis et aP (Il999i) . From this the CO column density is calcu- 
lated under the assumption that emission is optically thin and 
that the gas temperature is 100 K. The calculation of the col- 
umn density is insensitive to the choice of temperature, as long 
as it is in the interval of ~ 50-150 K. The calculation is more 
sensitive to the choice of velocity interval over which CO emis- 
sion is integrated, as this may vary by up to a factor of 2, when 
changing the velocity interval with 1 km s"' . The column density 
is then converted to total gas column density using a standard 
abundance ratio of C0:H2 of 10 '':1 and results are compared 
to CH3OH column densities to obtain an abundance. The inte- 
grated CO intensities, column densities and CH3OH abundances 
are all tabulated in Table |5] The gas column density in all three 
outflow positions is ~ 10^'' cm"^ resulting in abundances of ~ 
1-10x10"^ with respect to H2. Due to the uncertainties on the 
derived CH3OH column densities and on the gas column density, 
the abundances should be seen as order of magnitude estimates. 
Compared to the envelope abundances of SMMl and SMM4, 
this translates to methanol enhancements in the outflows of > 
10^. This is similar to what is found in other outflows; for exam- 
ple, th e enhancement for the LI 157 blue lobe is of the order o f 
~ 400 ('Bachi ller et al.lll993l : lBachiller & Perez Gutierrezlll997h . 
It should be noted, however, that the methanol line profiles are 
highly asymemtric in the LI 157 outflow as are the CO line pro- 
files. 

C'^O is often used as a tracer of quiescent gas column den- 
sity. If C'**0 J-3-2 emission is used rather than '^CO emission 
as above, the gas column density is increased by ~ two orders of 
magnitude to ~ 10^^ cm"^. This leads to a decrease in CH3OH 
abundance and therefore also a decrease in the enhancement. 
Representing the other extreme, typical enhancements are of the 
order of 5-50 (see Table |5]). Ideally, emission should be inte- 
grated over the same velocity interval for comparison, but in this 
case, part of the CH3OH line overlaps with the '^CO profile, 
part of it overlaps with the C'^O profile, which is why these two 
abundance measurements should be seen as extreme values. 

4. Discussion 

The inferred abundances are all larger than what can be produced 
by pure gas-phase reactions (e.g. Garrod et al. 2006). Thus the 
observed CH3OH gas must be formed on grains and subse- 
quently have desorbed. In the following, possible scenarios for 
methanol desorption and excitation are discussed. Finally spec- 
ulations on the nature of the individual young stellar objects in 
Serpens are made. 

4.1. CI-I3OI-I grain surface formation 

Methanol is formed on the surfaces of interstellar dust grains 
through hydrogenation of CO, a mechanism which has been 
studied in detail in t he laboratory and in theoretical models 
of surface chemistry ( Hiraoka et al J 120021 IWatanabe & Kouchil 
[2002; Fuchs et al. 2009). To form large amounts of CH3OH, it is 
imperative that CO freezes effectively out onto (water-ice cov- 
ered) dust grains. Toward several pre-stellar c ores (e.g., B68, 
Bergin et al. 2002) and Class objects (e.g., |j0rgensen et al.l 
2005 c) CO has been observed to freeze out very efliciently at 
temperatures lower than ~ 20 K and densities greater than ~ 10^' 
cm This "catastrophic" CO freeze-out has been observed di- 
rectly through observations of CO ice abundances which show 



an increase in the a mount of solid CO w ith respect to H2O ice in 
the densest regions (lPontoppidanl l2006) . Through use of Dusty 
modelling discussed above it is possible to estimate the extent 
of the CO freeze-out zone in each of the Serpens Class objects 
quantitatively, which is illustrated in Fig. |9] Here the predicted 
density and temperature profiles are shown as a function of dis- 
tance from the protostar itself. The density profiles of all four 
envelopes are very similar, however SMMl stands out in terms 
of temperature profile. Due to the higher luminosity (30 Lq) of 
SMMl with respect to the other source luminosities (~ 5 Lq), 
the envelope temperature is also higher throughout. 

The zone over which CO freezes out in SMMl is signifi- 
cantly smaller than that of the other sources, SMM3, SMM4 and 
S68N. In the case of SMMl, the present freeze-out zone extends 
from ~ 1000-5000 AU whereas the other envelopes have freeze- 
out zones starting at a few hundred AU and extend out to the 
same distance as for SMMl (see Fig.|9]i- This indicates that at 
present CO is not freezing out efficiently in the SMMl enve- 
lope, so any efficient methanol formation has probably ceased 
at this point in time. The colder envelopes surrounding SMM3, 
SMM4 and S68N could still be forming CH3OH. 

In terms of observi ng methanol directly in the ice itself, 
iPonto ppidan et al.l (l2004i) studied a small region extending south 
of SMM4 from 4000 AU to 12 000 AU by observing the 3.54 yum 
CH3OH features in absorption against background and embed- 
ded stars. The extent is illustrated by the white line in Fig. [T] 
In this region the CH30H-ice abundance is constant at 28% 
with respect to water ice or 3x10"^ with respect to gas-phase 
H2. This corroborates the interpretation presented here, that the 
gas phase abundance of CH3OH is low and constant out to 
12 000 AU in the SMM4 envelope. Beyond this line the CH3OH- 
ice abundance dro ps by at least an order of magnitude and 
'Ponto ppidan et al.l could only determine upper limits. It is in- 
teresting to note that the location where the ice abundance drops 
is where one of the outflow knots starts (SMM4-S; see Fig.[T]i. 
Thus the reason for the drop is a combination of the envelope be- 
ing more tenuous (riu < 10^ cm"^) far from the protostar, so that 
CO does not freeze out very efficiently, while at the same time 
whatever methanol is in the ice is sputtered into the gas phase by 
the outflow. 

ICuppen et al.l (l2009h have recently studied the formation of 
CH3OH on ice surfaces using a Monte Carlo method, in which 
the gas-grain chemistry based on laboratory data is simulated 
microscopically over long time-scales. The limiting factors in 
producing methanol is the availability of both CO and H on the 
grain surfaces, and so results show that CO hydrogenates effi- 
ciently to form CH3OH, especially at temperatures lower than 12 
K where atomic hydrogen can be retained efficiently. After 10^ 
years CH3OH may form up to 100 individual mono-layers on 
the grain, comparable to that of water ice. Hence, the ice abun- 
dance of methanol in the outer parts of the envelope will depend 
strongly on the temperature. Because the envelopes of SMM3, 
SMM4 and S68N are significantly colder than that of SMMl, 
they would still be actively producing methanol, which could 
explain the higher gas phase abundances (see Table|4]l. 

Observations show a different behaviour of C'**0 at the 
center position of S68N compared with the center position of 
SMMl. Figure [TOlpresents C'^O J-3-2 emission in colour and 
methanol emission from the 7o-6() A^ line overlaid as contours. 
The morphology of the C'^O emission shows a peak toward 
SMMl, as expected for a warm envelope, and a ring of emis- 
sion toward S68N. The C'^O abundances towards the two cen- 
tral positions are measured to be 3 and 1x10"^ respectively, 
based on integrated C'**0 intensities of 4. 3 and2.4Kkms"' and 
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Table 5. Gas column density, A^(H2) and CH3OH abundance in outflow knots as determined from '^CO 7=3-2 and C'^O 7=3-2 
emission. 
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Figure 9. Density and temperature profiles of the SMMl and SMM4 envelopes as obtained from Dusty modelling. The radial 
density profile is displayed in blue (left axis) and kinetic dust temperature in red (right axis). The profiles are shown for the SMMl 
envelope (full line) and the SMM4 envelope (dashed line). The CO freeze-out zone, characterised by temperatures lower than 25 K 
and densities greater than 10^ cm-^, is shown in dark gray, whereas the methanol evaporation zone (T > 80 K) is shown in light 
gray. To the right is a cartoon illustrating the differences between the envelopes surrounding SMMl and SMM4. Both illustrations 
are to scale and the dark blue regions correspond to the CO freeze-out zones while red indicates CH3OH evaporation zones. Note 
the much smaller CO freeze-out zone for SMMl. 



the gas column densities derived from dust continuum emission 
(Table|4|. This is to be compared to a standard C'^0 abundance 
of 1.8x1 Q-^ for a normal abundance of CO/ H? = IQ-"^ and an 
I6Q/I8Q (Wilson & Rood| [T994l) . Thus, CO is de- 

pleted by a factor of ~ 6 and 18 for the two sources as averaged 
over the entire envelope. For these two sources there is a clear 
anti-correlation between CO and CH3OH gas phase abundances. 

To examine whether this anti-correlation is unique to 
Serpens or a general feature of embedded source s, the outer- 
envelo pe CH3OH abundances of the sample of iMaret et al] 
(l2005h and lJ0rgensen et al.1 (l2005bl) was coupled wi fl i CO abun- 
dances for the same sources from 'j0rgensen et alJ (l2002h (see 
Fig- ED- The CO gas phase abundance averaged over the extent 
of the envelope is taken as a tracer of CO depletion in the sense 
that the total gas and ice CO abundance is assumed to be con- 
stant. For sources where the presence of a jump zone can repro- 
duce the observed CH3OH emission, only the outer abundance 
is plotted here. A typical uncertainty of 30% is assumed both 
for the CO and CH3OH abundances. Except for three sources 
(L1157, L1448-I2 and L483) there is a correlation between the 



two abundances as illustrated by the value of the Pearson corre- 
lation coefficient of r = 0.70. The Serpens sources are character- 
ized by the combination of low CO-abundance and high CH3OH 
abundance. This result implies that the CH3OH gas phase abun- 
dance is directly related to the current production of CH3OH 
in the outer, cold parts of the envelope, and that any difference 
from source to source is due to a difference in the amount of CO 
frozen out onto the grains. Thus, not all solid CH3OH is formed 
during the cold pre-stellar core phase, consistent with the lack of 
detected CH3OH ice toward background stars behind quiescent 
dense clouds. This also implies that the lack of CH3OH emis- 
sion at normal cloud positions between the sources is because 
the density is lower, and hence the timescale for CO to freeze out 
is much higher, i.e., the absolute CH3OH abundance will drop. 

4.2. CH3 OH desorption mechanism 

Once methanol has formed on a grain surface it can desorb ac- 
cording to two different mechanisms, thermal and non-thermal 
desorption. In the first mechanism the entire grain is heated ther- 
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Figure 10. Integrated C'^O, J -3-2 emission ( J Tmbdy) is shown 
in colour, with integrated CH3OH 7o-6() emission overlaid as 
contours. The image is centered on Serpens SMMl with S68N 
located to the north. Outflow directions are indicated by dashed 
(red-shifted) and dotted (blue-shifted) lines. 



mally (macroscopic grain heating), and the icy mantle evapo- 
rates entirely, releasing all adsorbed species into the gas phase. 
The ice mantles typically evaporate at grain temperatures of 
30-100 K, depending on species. The other mechanism is non- 
thermal desorption, in which the ice mantle is "heated" on a 
microscopic (local) scale, eithe r due to absorption of a single 
UV-photon (lOberg et al.ll2009ch. impact of a c osmic ray particle 
(Leger e t al.lll985l:lHasepawa & Herbs3ll993l:lHerbst & Cupped 
[2OO6), or the bindi ng energy being rel eased from the formation 
of a new mole cule (Garrod et al.ll2007h or sputtering in outflows 
jJimenez-Serra et al.i.2008i) . 



4.2.1. Thermal desorption 

In the case of methanol, the thermal evaporatio n temperature has 
been determined experimentally to 80-100 K (Bro wn & Bolinal 
[2007; Green et al. 2009) as is also indicated in Fig.|9] The lower 
temperature corresponds to evaporation of a pure CH3OH ice 
and the higher to a mix of CH3OH and H2O. Because CO freezes 
out on top of a water ice, CH3OH is expected to be present in the 
same layer Thus the desorption temperature should be close to 
that of a pure CH3OH ice (~ 85 K), which is slightly lower, but 
not much, than that of CH3OH mix ed with H20- T his is also 
found to be the case experimentally (lBisschopll2007h . Thus, for 
the thermal-desorption mechanism to be active, it is necessary 
that all grains are heated to greater than 80 K, i.e., close to the 
protostar itself. From the physical structure models above it is 
predicted that the radius at which Tdust > 80 K is of the order 
of 50-100 AU or a.'2-0'.'4 for t/=230 pc, depending on source. 
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Figure 11. CH3OH gas abundance versus CO gas abundance for 
a sel ection of Class and I sources based on d ata and analysis 
from Sch oier et all (|2002|) : iMaret et alj (l2005h : lJ0rgensen etal] 
(2005b) and this work. The Serpens sources are in red and 
IRAS 16293-2422 and NGC1333-IRAS2A are in blue. Sources 
for which only upper methanol limits exist are marked by ar- 
rows. The two dashed lines indicate the Icr confidence limits of 
the power-law fit. 



leading to a beam-dilution factor of 1500-5000. Methanol emis- 
sion originating from very close to the protostar itself has pre- 
viously been observed for two low-mass sources, IRAS 16293- 
2422 (c/=125 pc; van Dishoeck et al. 1995; Schoier et al. 2003) 
and NGC1333 IRAS2A (t/=250 pc; IMaret etaLH 2005). Even 
though the regions have not been spatially resolved with single- 
dish data, it has been possible to infer the existence of "hot 
cores" based on observations of high-ZT lines {K > 3) and mea- 
sured rotational temperatures of > 80 K. This indicates that 
the emission arises in a compact, warm and dense region, con- 
sistent with it originating close to the protostar itself and also 
consistent with high-spatial resolution interferometer data (e.g., 
|j0rgenseli et a l. 2005a, 2007). For the sources presented here, it 
has only been possible to provide an upper limit on the inner 
abundance, which in most cases is a few xlO"^ with respect to 
H2. One notable exception is SMMl, where the upper limit on 
the inner abundance is close to the outer abundance, a result that 
will be discussed further below (Sect. l43] ). 



4.2.2. Non-thermal desorption 

Since the temperature in the outer envelope is significantly lower 
than 100 K, the primary desorption mechanism must be non- 
thermal. This is supported by the fact that the abundance only 
follows the column density, and does not appear to depend on 
temperature. No attempt is made to distinguish between the 
above mentioned mechanisms here. If desorption is induced 
by secondary UV photons from cosmic ray ionization of H2 
the gas phase abundance is expected to be lO^'^-lO"-' times 
the ice abundance (Oberg et al. 2009ab). In the case of SMM4 
where the CH30H-ice ab undance is ~ 3x10"^ with respect to H2 
(iPontoppidan et al]l2004 ) the UV-induced desorption alone can 
account for a gas phase abundance of lO^^-lO"**. This is remark- 
ably close to the observed gas phase abundance of ~ 2x10"**. 

Garrodet al. (2007) modelled the release of CH3OH from 
the grain surface into the gas phase by examining whether the 
release of binding energy would be enough to evaporate the 
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molecule into the gas phase. They found that the fraction was 
probably in the range of 1-10%, but could not pin it down any 
further. In the case of SMM4, results presented here combined 
with ice observations indicate that much less than 1-10% is 
desorbe d, since the gas/ice abunda nce is ~ 10 Results pre- 
sented in lHasegawa & HerbstI (1 19931) show that direct cosmic ray 
desorption is not an efficient desorption mechanism for tightly 
bound species like H2O and CH3OH when compared to the 
accretion timescale. In fact, they find that the difference be- 
tween the two rates is of the order of 10"^. Previous studies show 
that direct cosmic ray ind uced desorption is most efficient for 
volatile species (see, e.g., S hen et al1 l2004; Ro berts et al. I l2007h 
and for species adsorbe d on very small grains (r < 0.05 fim; 
IHerbst & Cuppenll2006l) . 



4.2.3. Shock-induced desorption 

In outflows the main desorption mechanism is sputtering of the 
grain mantle, i.e., impacts of high-temperature gas molecules 
and atoms, primarily H2 and He. This can in principle be com- 
pared to the direct cosmic ray desorption mechanism, except 
that the flux of H2 and He is orders of magnitude higher and 
the impact energy is much lower. The efficiency of this desorp- 
tion mechanism is clearly seen from the fifth column of Table 
|5] where the CH3OH abundance in the outflow positions is 
compared to the ambient envelope abundance. In particular, the 
SMM4-S methanol outflow abundance is ~ lO"*" with respect to 
H2, i.e., only a factor of 3 lower than the measured CHsOH-ice 
abundance in a region located ~ 15" (4 000 AU) to the north 
(see Fig. [Tl and Pontoppid an et al.l 12004"). If the envelope ice- 
abundance in this region were also 3x10"^ prior to the impact 
of the shock, this would indicate that one third or more of the 
CH3OH ice is sputtered from the grain mantle. However, the 
measured gas-phase abundance as determined from '^CO emis- 
sion is an upper limit as briefly discussed above. The true value 
of the abundance probably lies in the range of 10 4o 10 where 
the lower limit is obtained from C'**0 emission. 



4.3. The nature of the embedded YSOs in Serpens 

Two low-mass YSOs have been confirmed observationally to 
have the same chemical characteristics as high-mass hot cores, 
IRAS 16293-2422 and NGC1333-IRAS2A. These two sources 
are at the same evolutionary stage as SMMl based on classi- 
cal indicators, e.g., Tboi. The characteristics are (1) a large num- 
ber of detected, saturated, complex organic molecules and (2) a 
high rotational temperature for methanol (> 80 K). At the same 
time, neither of these two envelopes s how s igns of extended 
methanol emission (van Dishoeck et an il995L Kristensen et al. 
in prep.). IRAS 16293-2422 is a close binary system surrounded 
by a circum-binary envelope. The inner parts of this envelope 
are passively heated to temperatures > 80 K, which, along with 
small-scale shocks impinging on the inner wall of the envelope 
(Chandler et al. 2005), give rise to the observed hot-core signa- 
ture. NGC1333-IRAS2 is also a binary (possibly triple) system, 
but the projected separation is greater (~ 30"). It is believed that 
the same mechanisms are at play in NGC133 3 -IRAS2A caus- 
ing the hot-core signatures (iMaret etal J 120051 : lj0rgensen et al.l 
l2005bl) . In the following the Serpens sources, and in particular 
SMMl, will be compared to these two low-mass hot cores, and 
differences will be discussed. 



There are several possible explanations for the lack of "hot 
core" characteristics in the Serpens sources which may be cate- 
gorized in the following manner: 

1. Physical: No gas is present close to the source, or no hot gas 
is present. 

2. Chemical: Hot gas is present, but the methanol abundance is 
very low. 

3. Observational: Warm methanol is present close to the source, 
but either the extent is very small, or emission is optically 
thick. 

IChoil (l2009h recently reported that SMMl is a close binary 
with a projected separation of ~ 500 AU. The primary would be 
the sub-mm source associated with SMMl (SMMla), while the 
binary (SMMlb) would be less embedded. SMMlb can be asso- 
ciated with a source observed by Spitzer at wavelengths shorter 
than 24 //m, where SMMla is not detected. If the projected dis- 
tance is similar to the actual separation between the two sources, 
then this could have cleared the inner part of the envelope of gas, 
explaining the absence of hot gas close to the (sub-mm) source. 
However, both IRAS 16293-2422 and NGC1333-IRAS2A have 
been show n to have inner holes or cavities similar in size to that 
of SMMl dSchoier et al.l2OO2Hj0rgensen et al.l2005d) . therefore 
this cannot be the entire explanation. 

Recent millimeter interferometry of SMMl indicates that 
the region c lose to the protostar has been cleare d of gas. Both 
lEnoch et all (|2009) and ivan Kempen et all (|2009|) report the de- 
tection of a resolved disk surrounding SMMl . The disk is unusu- 
ally massive (~ 1 Mr) and has a modelled radius of up to 300 
AU (lEnoch etal .1120091) . while the inner 500 AU of the envelope 
have been cleared of gas. In such a dense disk, only the upper- 
most layers will be warm or hot, depending on the distance to the 
protostar, as the dust extinction grows rapidly towards the mid- 
plane of the disk. Thus, the column density of hot met hanol will 
be ver y low. Through SED modeling of the disk, lEnoch et al.l 
(l2009l) estimate that the inclination is 15°, i.e., nearly face-on. 
This is in conflict with observations of the 3.6-cm radio jet, 
which indicate that the outflow is moving very close to the 
plane of the sky (Rodriguez et al. 1989; Moscadefli et al. 200^. 
Moreover, the '^CO line profiles are very symmetric along the 
large-scale outflow, i.e., the molecular outflows are also moving 
close to the plane of the sky. This would indicate that the disk 
is seen close to edge-on. If so, then the beam-dilution will be 
very high implying that any hot part will not be detectable in our 
single-dish beam. The disk surrounding NGC1 333-IRAS2A is 
less massive compared to SMMl (~ 0.056 M^: |j0rgensen et aH 
2009) and is viewed closer to face-on resulting in less extinction 
and lower beam dilution. This implies that the high-ZT emission 
observed in IRAS 16293-2422 and NGC1333-IRAS2A is origi- 
nating in the disk close to the protostar. The actual heating mech- 
anism (passive heating or small-scale shocks) cannot be distin- 
guished with the current observations. 

If the inner part of the envelope has not been cleared out 
of warm gas or if the disk is seen face-on, then it is possible 
that the abundance of methanol is low. As has been shown by 
current observations, the CH3OH gas-phase abundance is com- 
paratively high in the outer parts of the envelopes, so to de- 
crease the abundance in the inner part of the system, destruc- 
tion of CH3OH must be present. This destruction mechanism 
must be very efficient if it is to destroy all CH3OH in the hot- 
core parts of the envelope, where the abundance is expected to 
rise to > 10"^ w.r.t. H2. CH3OH can be destroyed in the gas 
phase through direct reactions with other species, however the 
rate coefficients are typically low. An alternative mechanism of 
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both the gas and the ice destruction is UV-dissociation. However, 
this destruction mechanism must also be at play in IRAS 16293- 
2422 and NGC1333-IRAS2A and there is currently no reason 
why UV-photodissociation would be more efficient in Serpens 
sources than in the other two nor that the UV-field is enhanced 
inSMMl. 

Finally, a low-mass hot core may be present in all of the 
Serpens sources, but not observable. This can be due to beam 
dilution or due to optical depth effects. The region over which 
methanol desorbs from the grain surface is expected to be ~ 
100 AU at most, and so beam dilution would be of the order 
of 10^. However, the same beam dilution would apply to at least 
NGC1333-IRAS2A, which is located at a similar- distance of 250 
pc, and can thus not be used as an argument. The line optical 
depth has been calculated in the Ratran simulations, and is typi- 
cally of the order of 0. 1 or less for the transitions observed here, 
even for lines arising in the inner-most part of the envelope. The 
dust opacity at 338 GHz is less than 0.08 at all times, as esti- 
mated from the Dusty modelling. 

Of the three explanations presented above the first is the 
more plausible if the distance is indeed 230 pc as assumed. 
Recent VLBA observations indi cate that the distance may be 
closer to 415 pc (iDzib et al.l2010 ). in which case the third expla- 
nation is more plausible. The other two reasons can be disproved 
through comparison with IRAS 16293-2422 and NGC1333- 
IRAS2A. However, only SMMl has been suggested to have a 
massive disk, and it is also the source that resembles the two low- 
mass hot cores the most in terms of luminosity. The other sources 
(SMM3, SMM4 and S68N) have considerably lower luminosi- 
ties by factors 4-10. Thus it may be that the hot-core regions 
around these sources are indeed much smaller and beam-diluted 
(~ 5x 10^), c omparable to sever al of the low-mass sources i n 
the sample of iMaret etaTI ( l2005h and lJ0rpensen et al.l ( l2005bl) . 
With the current observations only upper limits have been deter- 
mined of the molecular abundance in the inner envelope of the 
lower-luminosity sources. 

5. Summary 

Maps of rotationally excited methanol in the Serpens Molecular 
Core have been presented. Emission arises from the molecular 
envelopes of four deeply embedded sources and their associ- 
ated outflows. In particular, three outflow knots have been iden- 
tified based on their strong methanol emission. Mapping shows 
methanol emission in all of the four envelopes to be extended 
out to ranges of ~ 10000 AU. The abundance is constant in 
the outer parts of the envelope with a value of 10"''-10 de- 
pending on source. The methanol abundance at outflow positions 
is enhanced by up to 2-3 orders of magnitude with respect to 
the ambient abundance. The measured envelope abundances are 
consistent with non-thermal desorption of solid-state methanol, 
through, for example, UV-photodesorption. At outflow positions 
the enhanced abundance may be explained by sputtering of the 
grain mantles. The symmetric, slightly shifted line profiles point 
to the outflows being either caused by J-type shocks or moving 
along the plane of the sky. 

The CO gas abundance is found to be anti-correlated with 
the methanol gas abundance in the Serpens maps. This result has 
been extended to literature data of a large sample of other, simi- 
larly deeply embedded sources. The reason is that the more CO 
is frozen out from the gas phase, the higher the CO ice abun- 
dance, and this adsorbed CO is then converted to CH3OH ice 
on the grain surface from where it desorbs non-thermally. The 
non-thermal desorption mechanism implies that the gas-phase 



abundance follows the solid-state abundance closely, something 
which has directly shown to be the case here. Thus, the differ- 
ences in abundance between the four sources can be directly re- 
lated to the current production of CH3OH and reflects how much 
CO is frozen out now rather than in some possible colder past. 

The Serpens sources do not contain the chemical signatures 
of low-mass hot cores such as IRAS 16293-2422 and NGC1333- 
IRAS2A. In the case of SMMl, the more luminous of the 
Serpens sources and resembling most closely the two low-mass 
hot cores, the most likely explanation is that it harbours a mas- 
sive disk which is seen close to edge-on. The other three Serpens 
sources are all lower in mass and luminosity, and here the ab- 
sence of a hot-core signature could be ascribed to beam dilution, 
as also holds for the majority of the deeply embedded sources. 

In conclusion, the CH3OH emission is found to trace the fol- 
lowing: (1) energetic input into cold gas, primarily through out- 
flow interaction; (2) column density of cold gas in the outer en- 
velope; (3) reactions and subsequent desorption of grain surface 
products. 

Because water closely resembles methanol in the sense that 
it is exclusively formed on grains and it desorbs from grain sur- 
faces at r > 100 K, it is expected that the abundance structure 
will resemble that of methanol. This is to be tested with upcom- 
ing water-observations of these sources as part of the Herschel 
Guaranteed Time Key Project "Water in Star-forming Regions 
with Herschel" (WISH). 
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Table A.l. SED points used for Dusty modelling. 
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